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Abstract 

A Monte Carlo model to simulate nuclear collisions in the energy range going from SPS to 
LHC, is presented. The model includes in its initial stage both soft and semihard components, 
which lead to the formation of color strings. Collectivity is taken into account considering the 
possibility of strings in color representations higher than triplet or antitriplet, by means of 
string fusion. String breaking leads to the production of secondaries. At this point, the model 
can be used as initial condition for further evolution by a transport model. In order to tune 
the parameters and see the results in nucleus-nucleus collisions, a naif model for rescattering 
of secondaries is introduced. Results of the model are compared with experimental data, and 
predictions for RHIC and LHC are shown. 
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1 Introduction 



With the announcement of the discovery of Quark Gluon Plasma (QGP) at the Super 
Proton Synchrotron (SPS) at CERN the experimental heavy ion program moves 
now to the higher energies of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
and the Large Hadron Collider (LHC) at CERN. Whether this claim can be considered 
conclusive or not (see e.g. the most compelling experimental findings at the SPS 
0, 1^, |, 0, H, 3 are interpreted as positive signatures of QGP only when conventional, 
non-QGP models fail to reproduce them. Therefore, even in the case that QGP has 
already been obtained, it is most important that conventional models employed at the 
SPS become generalized for RHIC and LHC: They can be used to describe collisions 
between less massive nuclei or more peripheral events than those in which QGP is 
expected, and to establish the background to events with QGP production. 

On the other hand, the situation with conventional models is not clear at all. The 
description of a high energy collision between heavy ions is a complex task which 
involves different physical aspects. Predictions from different models for RHIC and 



LHC are far from being compatible, see the reviews ||T0[ and |[TT|]. For example, the 
values for central rapidity densities of charged particles coming from different models 
lie in the ranges 600 -ir 1500 for central AuAu collisions at RHIC and 2000 -r- 8000 for 
central PbPb collisions at LHC. 

In this paper a non-QGP model for collisions between nucleons or nuclei in the 
energy range going from SPS energies (~ 20 GeV per nucleon in the center of mass) to 
LHC energies (5.5 TeV per nucleon in the center of mass) is presented (different steps in 
this direction can be found in |12, |T^). The model is based on the ideas of Dual Parton 
Model (DPM) [g or Quark-Gluon String Model (QGSM) considering both soft 
and semihard components on a partonic level. These elementary partonic collisions 
lead to the formation of color strings. Collectivity is taken into account considering 
the possibility of strings in color representations higher than triplet or antitriplet, by 
means of string fusion, as done in |jl2|, (see related approaches in [jl^, |T^). String 
breaking leads to the production of secondaries. In this form, the model can be used as 
initial condition for subsequent evolution using a transport model, as those of |17|. 
Nevertheless, in order to tune the parameters of the model and apply it to nucleus- 
nucleus collisions, rescattering between secondaries is considered on the basis of 2 — > 2 
collisions, using a very simple model which allows us just to estimate the effects of such 
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process. The results of the code turn out to agree reasonably with existing experimental 
data on total multiplicities, and longitudinal and transverse momentum distributions, 
and semiquantitatively with strangeness production and stopping power. 

The paper is organized as follows: In Section |^ string formation will be discussed, 
both for soft and semihard components, whose separation will be established. Also in 
this Section collectivity, considered as string interaction or fusion, will be presented. 
Hadronization of the produced strings will be formulated in Section ^. In Section ^ our 
simple approach to rescattering between secondaries will be presented. A comparison 
with experimental data will be done in Section ^, and predictions for RHIC and LHC 



shown in Section o, together with some discussion on the first RHIC data [19, 20 



In the last Section we will summarize our conclusions and briefly compare with other 
approaches. 



2 Initial stage 

2.1 Elementary partonic collisions 

To compute the number of elementary partonic collisions we have to generate the 
partonic wave functions of the colliding hadrons. The steps to generate this wave 
function for the projectile A and target B are the following: First, the impact parameter 
b of the collision is generated uniformly between and Ra + Rb (in the case of nucleons, 
the total cross section determines the corresponding radius). 

Second, the nuclear wave function is computed. Nucleon positions inside the nucleus 
are distributed in transverse space according to a Woods-Saxon distribution for A > 11, 

l + exp[(r-r„)/a] ^ 
with Tn = 1.07A^/^ fm and a = 0.545 fm, and according to a Gaussian distribution for 



A < 11, with parameters chosen for each nucleus Then, Fermi motion is given 

to the nucleons in the nuclei uniformly in the range < p < pp-, with the maximum 
Fermi momentum given in the local Thomas-Fermi approximation ||2^ by 



PF = h [37rp(r)]V3, (2) 

with h = 0.197 fm GeV/c. 
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Now partons are generated inside each nucleon. Its number is given by a Poisson 



distribution 23 



iy^ocexp(-,(.)C/2) ""^l^'^ljl^ 9{s) = g,s'^l\ 9o=^-^, (3) 

with A = 0.139 the pomeron intercept minus 1, C = 3.0 the quasieikonal parameter 
which takes into account low mass nucleon dissociation, jp = 1.77 GeV^ the pomeron- 
nucleon vertex, ap = 3.3 mb the parton-parton cross section and y/s the center of mass 
energy for each nucleon-nucleon collision. 

Parton positions in transverse space (inside a nucleon) are given by a Gaussian 
according to Regge theory, 

^^^^ = 4^ '""P ("ix)' ^ = ^0 + 
with i?Q = 3.18 GeV^ and a' = 0.21 GeV^ the pomeron slope. 

Now, one parton from the projectile and one from the target produce an inelastic 
collision if both are within an area in impact parameter equal to ap = 27crp, rp = 0.23 
fm. In this way, events with no inelastic collisions are elastic, while those with at 
least one inelastic collision are inelastic. Taking the total cross section given by the 
quasieikonal model p3| 



(^tot = (yp f{z/2), z = — - — , f{z) = ^ _ ^, , (5) 

all cross sections can be computed, see next Subsection (all formulae reduce to the 
usual eikonal ones with C = 2). 



2.2 Semihard component 

The inclusion of semihard components, in the form of a two-component model, is 
needed to reproduce the pr spectra in hadronic collisions, see Section |^. In the model 
this is performed considering that an inelastic collision is hard with probability 

with Ah = 0.50, ^/so = 25 GeV and Ch = 0.0035. A hard collision proceeds through 
the packages PYTHIA 5.5 + ARIADNE 4.02 + JETSET 7.3 ^ Only gluon- 
gluon collisions are included in PYTHIA, and the key parameter here is the cut-off in 
transverse momentum 

PT™n = 3.03 + 0.11 In (s/so) GeV/c. (7) 



The minimum energy for an elementary collision to be accepted by PYTHIA is 20 
GeV, and for the global collision the minimum center of mass energy per nucleon is 
= 25 GeV. An event is considered hard if at least one of its inelastic elementary 
collisions successfully proceeds through PYTHIA. 

While the concrete choice of the parameters in pTmm comes from a fit to experimen- 
tal data, let us make some comments on its functional form. In our case, an increase 
of pTmin with increasing energy makes possible a smooth transition from the soft to 
the semihard part of the pt spectrum. Usually PTmm is taken as either constant or 
increasing as a polynomial of a logarithm of s p7| , [28[] . It may be argued that the 
PTmin valuc which indicates the transition from nonperturbative to perturbative QCD 
(pQCD), is related with the proposed saturation scale Ql 0: below this Q^, the 
number of partons in the hadron wave function cannot grow, as new partons fuse with 
the existing ones and cannot be resolved individually. Nevertheless, apart from con- 
ceptual differences, the dependences of PTmin and are not the same: while the first 
depends only on energy, the second one also depends on the size of the colliding objects 
{Ql oc A", a = 1/3^2/3). 

Results of the model for the total, inelastic (production) and hard cross sections 
in pp and pp collisions at different energies are shown in Fig. || and compared with 



experimental data for the total cross section It can be observed that both the 

total and the production cross section are too small at low energies, while they get 
reasonable values at higher energies. The reasons for the existing discrepancies are 
three: In first place, diffraction is not properly included in the model, so it is difficult 
to distinguish between production and inelastic cross sections. In second place, no 
reggeon contribution (decreasing with energy) has been included. In third place, at 
the level of the cross sections no distinction is made between nucleons and antinucleons 
as projectiles and targets. These two last reasons should improve the agreement with 
data at energies of SPS and Intersecting Storage Rings (ISR). Also in this Figure it is 
shown the value of PTmm and the mean number of total and hard inelastic collisions 
per event. 

2.3 String formation and fusion 



Each soft parton-parton collision gives rise to two strings [0, |T5[, stretched either 



between valence quarks and diquarks (for the first collision suffered by a nucleon) 
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or sea quarks and antiquarks (for the subsequent ones). For the latter, their flavors 
follow the ratio u : d : s = 1 : 1 : 0.26. Hard collisions proceed through PYTHIA as 
gg — > gg. For the string ends and hard gluons, the longitudinal momentum fractions 
are distributed as 

n 

L{xl,X2,...,Xn) = fgqiXl)fg{x2)---fg{Xn), = I- (8) 

k=l 

For soft strings ends, the individual momentum distributions are those of the QGSM 



fgg{x)=X^/\ /g(,-)(x) = ^, (9) 



with a lower cut-off Xmin = 0.3 GeV/ ^/s^n to ensure that the strings have mass enough 
to be projected onto hadrons, ^/snn being the center of mass energy per nucleon. 

For partons involved in hard collisions, the longitudinal momentum fractions are 
taken by PYTHIA from PDFLIB |3^, with the possibility of considering the difference 
of parton distributions inside nuclei given by the parametrization EKS98 or by 



a parametrization as F2A [HI- After generating the final gluons, each of them splits 
into a {qq) pair and strings are stretched between them, according with the standard 



procedure in PYTHIA [|25|. 

The transverse momentum of both partons at the string ends and hard partons, 
coming from a nucleon which has been wounded m times, is given by a Gaussian: 

Tip±) = 4^ exp {-pI/6^), 6 = 0.5^ GeV/c; (10) 



in this way, pT-broadening is taken into account |35 



The number of strings exchanged in one collision is quite low for nucleon-nucleon 
collisions, but this number increases with energy, size of projectile and target and 
centrality of the collisions. Strings can be viewed as objects with a certain area, given 
by the uncertainty relation as oc 1/ (pf^) , in the transverse plane of the collision. When 
the number of strings is high enough, they begin to overlap and the usual hypothesis 
in QGSM or DPM of the strings being independent sources of secondary particles is 
expected to break down. A possible way of considering this is to compute the density of 
strings in the transverse plane and use two-dimensional percolation as an indicator of 
the onset of collectivity [p6| , p^ . Percolation takes place when domains of overlapping 
strings acquire a size of the order of the total available size for the collision. 

While percolation is a second order phase transition, the option we use in this 
model, fusion of strings, does not lead to any phase transition In the model. 
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ordinary strings (i.e. in a triplet representation of SU(3)) fuseQ in pairs when their 
parent partons (those which determined the inelastic collision the strings come from) 
are within a certain area afus = 27rrj„^ in impact parameter space. In the code we 
consider only fusion of two strings but there is a probability of fusion of more than 
two. An effective way of taking this into account is to increase the cross section for 
the fusion of two strings, for which we will take cxp < a jus = 7.5 mb {rf^s = 0.35 fm). 
This value is crucial to reproduce the strangeness enhancement in central SS and SAg 
collisions at SPS The fusion can take place only when the rapidity intervals of the 
strings overlap. It is formally described by allowing partons to interact several times, 
the number of interactions being the same both for projectile and target. 

The quantum numbers of the fused strings are determined by the interacting par- 
tons and their energy-momentum is the sum of the energy-momenta of their ancestor 
strings. The color charge of the resulting string ends is obtained according to the SU(3) 
composition laws: 

{3}® {3} = {6}©{3}, {3}®{3} = {1}©{8}. (11) 

Thus, two triplet strings fuse into either a sextet or an antitriplet string with probabil- 
ities 2/3 and 1/3 respectively, and one triplet and one antitriplet string fuse into either 
a singlet or an octet string with probabilities 1/9 and 8/9 respectively. 

Two comments are in order: On the one hand and as written above, the fusion 
of strings means nothing related to a phase transition. On the contrary, percolation 



of strings |36] is a non-thermal second order phase transition. In this case, the key 
parameter is 1] = vrr^ A^/ (tt/?^) , which is the density of strings N/{ttR\) (number of 
strings A^ produced in the overlapping area of the collision, ttR^ for central collisions) 
times the transverse size of one string vrr^. The critical point for percolation is t]c — 
1.12 -f- 1.5 depending on the profile function of the colliding nuclei With r ~ 



0.2 0.25 fm, this critical value means 6^12 strings/fm^. The value of 9 is reached in 
central PbPb collisions at SPS, in central AgAg collisions at RHIC and in central SS 
collisions at LHC We expect for rj around or greater than rjc, that the approximation 
of fusion of just two strings fails. On the other hand, only fusion of soft strings is 
considered. Hard strings are not fused, their area being proportional to l/pr- Some 
effect of the fusion of such strings could appear at LHC energies where, for instance, in 
central PbPb collisions they amount for 32 % of the binary nucleon-nucleon collisions, 
similar mechanism exists in RQMD jl^], called color ropes. 



3 Hadronization 



Now we consider the breaking of a soft string with color charges Q and Q in its ends 
(corresponding to a representation {A^} of SU(3)). In our model, it is due to the 
production of two (anti) quark complexes with the same color charges Q and Q as 
those at the ends of the string |T^0. The probability rate is given by the Schwinger 



formula 40 



W cc Kf^y exp (-vrMl/iTi^v}), 



(12) 



where K^^y is the string tension for the {A^} representation, proportional to the cor- 
responding quadratic Casimir operator C^j^iy (as found both in lattice QCD and in the 



Stochastic Vacuum Model 41, 42 
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(13) 



For the longitudinal breaking of the string, an invariant area law is employed 



P oc exp (-M), K{N} oc bCf^y, A = p+p 



(14) 



being the area in light-cone momentum space determined by the breaking point in the 
center of mass frame of the string. This law gives results quite similar to those of the 
Lund model |^ implemented in JETSET p5 |. 



We proceed as follows: Eq. ([T^ ) is used to decide the flavors of the quark and 
antiquark complexes created. We take K^^y = 0.18 GeV^ and rriu = rrid = 0.23 GeV/c^, 
rUs = 0.35 GeV/c^, and the masses of a complex (gi . . . qi) is given by M(gi . . .qi) = 
Yli=i f^gi- Then px is given to one of the created complexes and —px to the other one, 
according to a Gaussian law 



f{Px) oc exp(-a|;v}PT)) 



(15) 



with a 



{3} 



a 



{3} 



4 GeV-2 and 



a 



{TV} 



2a 



C 



{3} 



{3} 



C 



, {iV}^{3},{3}. 



(16) 



{TV} 



^This possibility is the dominant one for strings formed by fusion of two triplet strings [|T2|. For 
higher color representations, production of quark complexes with color charges Q' < Q begins to 
dominate. This option is taken in RQMD nevertheless, the close similarities in the consequences 
of string fusion in both approaches, strongly suggest that the difference in the breaking mechanism 
can be compensated by a different choice of the fragmentation parameters. 
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Finally a breaking point is sampled according to Eq. (|T4D in the available phase space, 
with b = 1.83 GeV'\ 

Fragmentation proceeds in an iterative way: String fragments are taken as new 
strings which are broken again, until the mass of the created fragments is too low 
to allow further breaking (i.e. projection onto hadrons with the right quantum num- 
bers). Then these final fragments (and those fused strings resulting in the singlet {1} 
representation) are treated as quark clusters and decayed according to combinatorics 
and phase space. Spin of the produced particles is constructed according to SU(2) 
considerations. 

The main consequences of string fusion are a strong reduction of multiplicities (both 
due to the energy-momentum conservation and to the reduction of the effective number 
of sources of secondaries) and a slight increase of (p^) [|T^, an increase in baryon and 
strangeness production [jT2|, a strong increase in the cumulative effect and a 



decrease in forward-backward correlations \i6 



On the other hand, strings produced in hard collisions (only gg — > gg) are man- 
aged by PYTHIA + ARIADNE + JETSET |2§g. For ARIADNE, PARA (6) is fixed 
so as the transverse momentum of the radiated gluon should be less than that of the 
hard gluon (i.e. the one participating in the gg scattering) and MSTA(9)=MSTA(14) = 
MSTA(31)=0. In JETSET, PARJ(41)=1.7 GeV-^ and PARJ(42)=0.6 GeV"^. besides, 
PARJ(21)=0.55 GeV/c. Only production of three flavors (u,d,s) is considered in the 
present implementation of the model, so it cannot be used to study production of heav- 



ier flavors (some steps in this direction were done in [|T3|). As a last point, MSTU(4), 
MSTU(5) and the dimensions in LUJETS have been set to 120000, which has shown 
to be enough for central PbPb collisions at LHC All the other parameters and options 
in these programs have been set to their default values. 



4 Rescattering of secondaries 



As stated in the Introduction, in this stage the model can be used as an initial con- 
dition for further evolution, using either a hydrodynamical model or a microscopic 

. (see 



transport as RQMD pi, UrQMD [Tl, HSD ART p 



13] for a study of 



evolution of particle and energy densities). Nevertheless, it is usually assumed that the 
■^Although it is not relevant for the generation of momenta of final gluons, in PYTHIA and JETSET 
the value of Aqcd has been taken from the corresponding set of parton densities for 5 flavors. 
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enhancement of hyperons, antihyperons and 0's observed in heavy ion experiments at 
SPS [|, H, H, 0, H cannot be fully explained by using exclusively a mechanism which 
goes beyond the independent string hypothesis, as string fusion [IB], |T7|, ^ or 
baryon junction migration [18, ^ |9, Q In order to reproduce these experimental 
features rescattering of particles in the hadron gas (produced particles among them- 
selves and with spectator nucleons) has been introduced in many models. To tune 
the code and study nucleus-nucleus collisions, we will make a very simple rescattering 
model with no space-time evolution, fitted to SPS data. Results of this approach will 
be presented, but one must keep in mind that predictions which depend critically on 
rescattering effects should be taken with much care. 

Our implementation of rescattering is extremely naif, trying not to solve the full 
Boltzmann transport equation for all particles but only to make a model, as simple 
as possible, which gives us an estimation of which effects such rescattering would pro- 
duce. Neither formation time nor space-time evolution of secondaries are properly 
considered; instead we require a common minimum density of particles in the rapid- 
ity bin of the considered particles, for rescattering to occur. This minimum density, 
dN / {dydpT)\min = 17, has been chosen for rescattering not to affect results in nucleon- 
nucleon collisions up to the highest energies. Rapidity and pt distances between par- 
ticles have to be lower than 1.5 units and 0.3 GeV/c respectively. Only two body 
reactions have been included, with inverse reactions as required by detailed balance. 
Spin is ignored, and rescattering takes place before resonance decay. All cross sections 
are taken equal for all reactions (except for Q production and nucleon annihilation). 

Operationally, both products of string breaking and spectators are randomly or- 
dered into an array (1, . . . , N). We compute the possibility of rescattering of the first 
element with all the others in pairs: (1, 2), (1, 3), (1, 4),. . . If either in one of the pairs 
(l,j), j = 2,...,N, rescattering (either elastic or inelastic) occurs or we reach pair 
(1, A^) with no kind of scattering happening, we go to element 2 and examine the pairs 
(2, 1), (2, 3), (2, 4),. . . This is repeated until the pair (A^, A^ — 1) is examined. As parti- 
cles produced in rescattering of the pair (i, j), i,j = l,...,N,ij^j, occupy the same 
places i,j in the array as their ancestors, particles produced by rescattering have a 
chance to rescatter again. 

The probability for two particles to scatter in a given inelastic channel is 7 % 
(except for channels involving fi's and nucleon annihilation, where it is 14 and 70 % 
respectively). For a given process, the probability for elastic scattering is given by 
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the sum of those corresponding to all inelastic channels considered for these initial 
particles. Cross sections (probabilities) are considered energy independent, except for 
the trivial kinematical thresholds, and isotropic in the center of mass of the colliding 
secondaries and/or spectators. The considered reactions (together with those for the 
corresponding antiparticles) can be classified into: 



Light pair, (gg), annihilation to create another light pair, or light quark exchange: 
TiN ttN, 7T71 TTTT, TiY ttY , ttH ttS, KN KN and KY — > KY, where 
F = S,A. 



• Other considered reactions are: vrA^ —>■ KY, vrvr — > KK, ttY KS, vrS Kfl 
and KN (pY. These reactions can be classified into: 

1. Light pair, [qq), annihilation to create a (ss) pair. 

2. Reactions with baryon exchange (that is, with three lines in the t-channel). 

• Reactions with strangeness exchange: KN ttY, KY ttH, KS nQ, 
KY —>■ (f)N, KS —>■ (j)Y and KVt (pY. This type of processes can produce 
(anti)baryons with several strange (anti)quarks and are exothermic. 

• Nucleon-antinucleon annihilation into two pions: NN tttt. This type of re- 
action has a much larger cross section at low energies than reactions consider 
before; for this reason its probability has been chosen ten times larger than the 
others. This is also an effective way to take into account final states involving 
more than two pions. 

To simplify, particles produced in rescattering are always projected onto the lowest 
spin state. Decay of resonances proceed through the usual JETSET routines, with 
MSTJ(22)=2, and decay of vr^'s is forbidden. The results of our rescattering model 
on strangeness and baryon/antibaryon production can be summarized in three points: 
hyperon and enhancement, antinucleon annihilation and a slight increase of stopping 
power (kinematical effects of our rescattering model are very small, due to the applied 
cuts in rapidity and transverse momentum). Besides, a slight decrease of multiplicities 
appears, as we will see in the next Section. 
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5 Comparison with experimental data 



In order to show the quahty of the choice of the parameters, in this Section we will 
compare the results of the code with experimental data. We will also analyze the 
influence of the different physical mechanisms implemented in the model. From now 
on and unless otherwise stated, results of the code come from its default version with 
string fusion, rescattering (which do not affect results in nucleon-nucleon collisions, and 



in pA collisions at SPS energies), and GRV 94 LO [§3| parton densities with EKS98 
nuclear corrections. 



5.1 Hadron-hadron collisions 

Results of the model for the mean numbers of produced particles in minimum bias pp 
collisions at y/s = 19.4 and 27.5 GeV are shown in Tables ^ and ^ respectively, com- 
pared with experimental data. An overall agreement can be observed. Two comments 
are in order: On the one hand, the influence of fusion in nucleon-nucleon collisions is 
tiny, apart from a slight increase in antibaryons. On the other hand, the number of 
both A's and A's is overestimated in the model. This is due to the fact that in the 
model, threshold effects, important at these low energies, are treated very roughly (see 
further comments in the next Subsection). At higher energies the situation improves. 
For example, in pp collisions at ^/s = 200 GeV, the mean number of A + A in the 



model is 0.56, to compare with the experimental result 0.46 ± 0.12 ||56[ . 

In Fig. 0, rapidity and transverse momentum distributions of negative particles in 
minimum bias pp collisions at ^/s = 19.4 GeV are shown and compared with experi- 
mental data. The result is satisfactory. In Figs. ^ and ^ pseudorapidity and transverse 
momentum distributions of charged particles in pp collisions at ^/s = 200 and 1800 
GeV are compared with experimental data. The agreement is reasonable, although the 
multiplicity at 200 GeV seems to be slightly underestimated. The results of the model 
without semihard component are also shown in these Figs., and they do not describe 
the pt distributions, which justifies the inclusion of hard collisions. Besides, the results 



given for different sets of parton distributions, both old ||5^ and new ||61| , |6^] and lead- 
ing order or next-to-leading order, are very similar. This fact may look surprising from 
a pQCD point of view. The main reason is that the cross sections and the number of 
inelastic collisions in our model are determined by Eqs. (|^), (^, (^ and (|^), which are 
independent of the choice of partonic distributions in PYTHIA (this is not so in other 
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models, see e.g. [^). We also think that the quite high prmin we use in PYTHIA, Eq. 
(0), and the gluon radiation and fragmentation performed by ARIADNE and JET- 
SET, may have some influence on the fact that no difference is apparently seen in the 
transverse momentum distributions. 

In Fig. ^ the evolution of the mean transverse momentum of charged particles is 
studied in pp collisions at SppS, versus the center of mass energy and, for different 
particles, versus central charged multiplicity. The trend of data is reproduced and we 
find the agreement reasonable (this cannot be achieved without the hard component, 
as seen in this Figure). In Fig. ^ the topological cross section for charged particles 
in the central region is examined at different energies for pp collisions at SppS, and 
the agreement is also reasonable, considering that the model slightly underestimates 
multiplicities at 200 GeV but correctly reproduces those at 1.8 TeV, see Figs. |^ and^. 

5.2 Proton-nucleus and nucleus- nucleus collisions 

In Table ^ results of the model in pA collisions are compared with experimental data 
on negative multiplicities. An overall agreement is obtained. The reduction of multi- 
plicities due to string fusion can be observed. 

In Table |, mean numbers of produced particles are compared with experimental 
data, for central SS collisions at SPS energies. The agreement is reasonable. Only the 
number of both A's and A's in the model is significantly below the experimental data. 
The number of A's is increased by both string fusion and rescattering, while that of 
A's is mainly determined by only string fusion (see results in PbPb below). Anyhow, 
rescattering is seen to have little effect in SS. 

Let us now discuss PbPb collisions at SPS. In the last year a large excitement has 

arisen in the heavy ion physics community, related to the possibility of Quark Gluon 

Plasma (QGP) already been obtained at SPS energies [|l|. In particular several signals 

were mentioned, which point out to the existence of QGP. Putting aside the abnormal 

J/tp suppression and the excess of dileptons found, there are three signals related to 

baryon and strangeness production, namely the large enhancement of the (anti)hyperon 

yields (A, H, Q) in PbPb collisions compared to pPb, observed by the WA97 |Q and 

the NA49 0] Collaborations^; the linear increase of the inverse exponential slope of 

recent reanalysis Q of S data done by the NA49 Collaboration gives yields at midrapidity 
which are in much closer agreement to the WA97 Q results than the previous analysis of NA49 . 
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the rriT distributions ('temperature') in PbPb collisions with the mass of the observed 
particle, except for Q |70|; and the different behavior of the temperature between pp 
and AA collisions. These characteristics have been interpreted as the existence of an 
intrinsic freeze-out temperature and a collective hydrodynamical flow which is gradually 
developed: firstly, for SS collisions, and, in a more clear way, in PbPb collisions. In this 
Subsection we will examine some of these points using our model, together with other 
interesting aspects as production 0, |^, different particle ratios and stopping 
power 1?^ , 

In Fig. 1^ we show our results for Q, H and A yields for pPb, and central PbPb 
collisions at SPS with four different centralities, together with the experimental data. 
In order to disentangle the different processes contributing, in Fig. |^ it is shown the 
results of the code for central {b < 3.2 fm) PbPb collisions without string fusion and 
rescattering, with string fusion, and with string fusion and rescattering. A reasonable 
agreement with data for PbPb is obtained, only the f2's are a 40 % below the data 
and we have some excess of A and S"*", see next paragraph. Similar results have been 
obtained in the Relativistic Quantum Molecular Dynamics model |T^, ^ by a mech- 
anism of color ropes which consider fusion of strings; also in the Ultra Relativistic 



Quantum Molecular Dynamics model |17| and in the HIJING model |Tq] by using an 



ad hoc multiplicative factor in the string tension. Also the Dual Parton Model |14 



considering the possibility of creation of diquark-antidiquark pairs in the nucleon sea, 
together with the inclusion of diagrams which take into account baryon junction mi- 
gration |4^, [5^, can reproduce the experimental data (for f2's some rescattering 
has still to be added). The string fusion is the main ingredient to obtain an enhance- 
ment of A production and also to reproduce the S data. However rescattering seems 
fundamental to get enough f2's. 

Nevertheless, our results for pPb are higher than the data for S"*" and A; this last 
feature looks quite strange, as A and S~ yields agree with data, but we overestimate 
both A and A production in pp collisions at this energy^, see Table |l[ As rescattering 
plays a minor role in minimum bias pPb collisions, this turns out to be a result of string 
fusion. About A, our results are higher than the WA97 data also in PbPb, its pro- 
^In our opinion, the comparison of (anti)hyperon nucleus-nucleus data with those in nucleon- 
nucleon collisions should be taken with caution at SPS, because at this relatively low energy the 
nucleon-nucleon value rises sharply with increasing energy due to the tmm- aud delayed threshold 
effects p3, which usually are not properly implemented in models. 
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duction being mainly determined by string fusion and hardly affected by rescattering. 
This fact makes that our results for PbPb are really an extrapolation in the model from 
the value for A production in central SS collisions by the NA35 Collaboration, which 
was used to fix the fusion cross section afus [^11 (even so, the model underestimates A 
production in central SS, see Table |^). So, from the point of view of our model, there 
exists either a large A annihilation or a conflict between NA35 data for SS and WA97 
data for PbPb and pPb. 

In Fig. ^ we plot the inverse exponential slopes of the distributions for different 
particles, together with the WA97 experimental dataQ. A semiquantitative agreement 
is obtained. In particular it can be seen that the Q slope does not obey the linear 
increase with increasing mass both in the model and in data, and that rescattering 
slightly increases temperatures. 

About (p enhancement, our integrated yields per event without fusion, with fusion, 
and with fusion and rescattering are 3.55, 4.20 and 5.35 respectively, in rough agreement 
with experimental data, 7.6 ± 1.1 @]. In Fig. ^ the stopping power is shown, i.e. 
the p — p rapidity distributions for central PbPb collisions at SPS, compared with 
the experimental data |72|, together with the predictions for RHIC and LHC energies. 
This quantity is essentially determined by the string fusion mechanism and rescattering 
only plays a minor role. As discussed for strangeness enhancement, it has been pointed 
out that baryon junction migration ^ |5ll will enhance the stopping power due to 
diagrams additional to the usual ones of the Dual Parton Model. The inclusion of these 
diagrams also explains the SPS data. We have not taken into account such diagrams 
to avoid double counting, because in the fusion of strings they are partially included 
in an effective way. In Fig. |ll] the antiproton rapidity distribution in central PbPb 



collisions is presented and compared to the experimental data ||7J]; a great suppression 
of the antiproton yield is seen, due to rescattering. 

In Table |] our results for the ratios between different particles are compared with 



the experimental data |]7T| for PbPb central collisions at SPS. We observe an overall, 
rough agreement with the SPS data, with some excess of A and S"*", see Fig. ^ and 
comments above. 

Let us emphasize that we obtain a semiquantitative agreement with the experimen- 



^The fits have been performed in the same ttit regions as WA97 did jT^l- For statistical reasons, 
we compare the slopes in the model for yields integrated over all rapidities, with experimental data 
taken in the central rapidity region. 
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tal data in PbPb, in three of the features advocated as signals of QGP production. 
We are only below data in Q production by less than a factor 2. So we think that our 
rescattering model, being very simple, can be useful as a tool to show the trend of such 
effect and at least help to tune the initial condition which can be used in transport 
models. 

Finally, let us comment on multiplicities in PbPb collisions at SPS energies. For a 
centrality of 5 % (corresponding in the model to 6 < 3.4 fm), we get, for dN~/dy at 
y = 0, 265, 250 and 235 without string fusion, with string fusion, and with string fusion 
and rescattering respectively. Experimentally, the NA49 Collaboration gets 196 ± 10 



75[| , while the WA97 Collaboration gets 178 ± 22 [Q. In view of these data the code 



overestimates multiplicities. On the other hand, if we compare the charged multiplicity 
per participant (wounded) nucleon and pseudorapidity unit at midrapidity versus the 
number of wounded nucleons in PbPb collisions at SPS, with data from the WA98 



Collaboration [^, the trend of data seems to be reproduced, while their magnitude 



is underestimated ||T8[. In Fig. |12| we show the rapidity distribution of negatives 
compared with NA49 data ||75|. 



6 Predictions for RHIC and LHC 

Predictions for pseudorapidity and transverse momentum distributions of charged par- 
ticles in nucleon- nucleon and central nucleus-nucleus collisions at RHIC {^/s = 200 
GeV per nucleon) and LHC {^/s = 5.5 and 14 TeV per nucleon for nucleus-nucleus and 
nucleon-nucleon collisions respectively) can be seenQ for nucleus-nucleus and nucleon- 
nucleon collisions respectively) in Figs. [1^, |1^ and |15[ While at SPS the influence of 
string fusion on multiplicities at midrapidity is of the order 10 -r- 15 %, at RHIC it 
reaches a 30 ^ 35 %. In these Figures, the large influence of the hard contribution at 
LHC can be observed. Again, the striking fact of the small influence of parton densities 
appears, both in nucleon and in nuclear collisions. On the other hand, the scaling of 
nucleon-nucleon with the number of wounded nucleons (the Wounded Nucleon Model 
Il79|| ) gives predictions which lie far below any of those of our model. 

In Fig. y we plot the inverse exponential slopes of the ttit distributions for different 
particles at RHIC. We see that, compared to the SPS situation, temperatures get higher 
^ About the reliability of predictions for RHIC and LHC, see comments in the last paragraph of 
Subsection 2.3 and in the first paragraph of Section ^ 
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in all cases, as expected. 

We present our predictions for different particle ratios at RHIC and LHC in Table 
It can be observed that our results are not very different to those of statistical models 



T0| , pQ , |81| , |8^ . However, strangeness enhancement in our case has nothing to do 
with thermal and/or chemical equilibrium. The main difference in the predictions for 
RHIC and LHC between the String Fusion Model and statistical models is the overall 
charged multiplicity, which is respectively 950 and 3100 for SFM and 1500 and 7600 for 
statistical models (assuming initial temperatures of 500 and 1000 MeV for RHIC 
and LHC respectively). Besides, predictions for the stopping power at RHIC and LHC 
energies are presented in Fig. |10|. Now, a pronounced dip appears at midrapidities. 

Detailed discussions on first RHIC results will be given elsewhere ||7^ . Here we sim- 
ply compare our results with some preliminary data of the PHOBOS [|T^ and PHENIX 
pm Collaborations at RHIC. For charged particles we obtain dN/drj ||r;|<i= 520 and 585 
for the 6 % more central AuAu collisions at ^/s = 56 and 130 GeV per nucleon respec- 
tively, to be compared with 408 ±12 (stat.)±30 (syst.) and 555 ±12 (stat.)±35 (syst.) 
(609 ± 1 (stat.) ± 37 (syst.)) in PHOBOS (PHENIX). Our prediction for = 200 
GeV per nucleon with the same centrality cut is dN/drj \ \r^\<i= 635. 

7 Conclusions 

A Monte Carlo model0 for nucleon and nuclear collisions in the energy range going from 
SPS to LHC has been presented. It is based on a partonic realization of Regge-Gribov 
and Glauber-Gribov models and its translation to strings following the DPM/QGSM 
ideas. A hard component is included to reproduce the high transverse momentum tail 
of the spectrum. Collectivity is included considering the possibility of fusion of pairs 
of strings. Strings are decayed in a conventional way. In order to tune the parameters 
of the model and apply it to collisions between nuclei, a naif model of rescattering has 
been introduced. 

The results of the models turn out to agree reasonably with total multiplicities, 
and longitudinal and transverse momentum spectra in the energy range from SPS to 



®The code, called psm-1.0, has been written in Fortran 77 and can be taken as 
an uuencoded file, containing instructions for installation and use, by anonymous ftp 
from [Ftp: / / ftp.uco.es / pub/falarpcn/ , or from the web sites http : / / www. uco . cs/ 'fal arpen / or 



http://fpaxpl.usc.es/phcnom/, or requested from the authors. 
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TeVatron. The agreement with strangeness production, temperature behavior and 
stopping at SPS is semiquantitative. 

There exist other Monte Carlo models for multiparticle production in nuclear col- 
lisions at ultrarelativistic energies (see [11| for a review): RQMD [16|, UrQMD |T7|, 



HIJING 13, DPMJET [|8l, HSD [|7|, NEXUS |8|], VNI g, AMPT |8|, LUCIFER 



8^,. . . Let us examine the main similarities and differences, concerning the stage be- 
fore rescattering is applied. Both DPMJET and our model are realizations of the 
DPM/QGSM which include a hard component, but we introduce string fusion, while 
DPMJET considers diquark breaking diagrams. RQMD takes into account string fu- 



sion (and now UrQMD and HIJING [0 in a simple way), but no hard part is included 
either in RQMD or in UrQMD. The main difference with HIJING lies in the soft 
component, which is considered energy independent is HIJING (and in this way, the 
multiplicity increase with increasing energy is mainly due to the hard component), 
while in our case it increases as an unitarized supercritical pomeron. VNI is a parton 
cascade code, in which the initial stage is mainly generated by hard collisions, with no 
hadronic degrees of freedom (strings). AMPT is a hybrid code, which uses HIJING 
as initial condition for parton cascade and, after hadronization, performs hadronic 
transport. HSD is focused in the transport of hadronic degrees of freedom, the ini- 
tial stage not coming from strings stretched between partons of projectile and target, 
but considering strings as excitations of nucleons in the projectile and target, as in 
Fritiof similar comments can be made for LUCIFER. Finally, NEXUS is based in 
DPM/QGSM, trying to solve the problem of energy-momentum conservation for both 
cross sections and multiparticle amplitudes at the same time. In our model, energy 
conservation is strictly taken into account only for multiparticle amplitudes. Besides 
NEXUS takes into account triple pomeron diagrams, which in our case are effectively 
included in string fusion. 

A detailed comparison of results of the model with the first RHIC results will be 
presented elsewhere |7^. As future developments, strangeness production should be 
reconsidered and production of heavier flavors included. Also fusion of more than two 
strings and the possibility of a phase transition like percolation of strings is needed in 
order to improve predictions for LHC and study the possibility of QGP formation in 
the framework of string models. 
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List of figures: 



1. Upper plot: Results in the model for the total (solid line), production (dashed 
line) and hard (dotted line) cross sections versus -y/s, compared with experimental 
data for total cross sections in pp (filled circles) and pp (open circles) collisions taken 
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mean number of total (dashed line) and hard (dotted line) inelastic collisions per event 
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contribution, which decreases quickly with increasing energy, is not included, see text. 

2. Results of the code for the rapidity distribution (upper plot), and the pr distribution 
for particles with 2 < yiab < 4 (lower plot), of negative particles in minimum bias pp 
collisions at piab = 200 GeV/c, compared with experimental data |5^, |58[ . 



3. Results of the code for the pseudorapidity distribution (upper plot), and the pr dis- 
tribution for particles with \ri\ < 2.5 (lower plot), of charged particles in minimum bias 
pp collisions at ^/s = 200 GeV, compared with experimental data [^, Solid lines 
are the results with GRV 94 LO parton densities [Q, dashed lines with CTEQ5L |^ 



dotted lines with GRV 98 HO ^2[ and dashed-dotted lines results without semihard 
contribution. In the pt distribution, model results have been normalized to experi- 
mental data. 

4. Results of the code for the pseudorapidity distribution (upper plot), and the pt 
distribution for particles with |?7| < 1 (lower plot), of charged particles in minimum 
bias pp collisions at -y/i = 1.8 TeV, compared with experimental data Line 



convention is the same as in Fig. ^. In the pt distribution, model results have been 
normalized to experimental data. 

5. Upper plot: results of the code without hard part (dotted line), without string fusion 
(dashed line) and with string fusion (solid line) for (px) of charged particles with |?7| < 
0.5 in pp collisions, versus ^/s, compared with UAl data and a parametrization 
given in this reference (dashed-dotted line). Lower plot: results of the code for (px) of 
tt"^ (solid line), K"^ (dashed line) and p (dotted line) in pp collisions at y/s = 1.8 TeV, 
versus the pseudorapidity density of charged particles for \ri\ < 3.15, compared with 
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E735 data ||66| for vr^ (circles), K"^ (squares) and p (triangles). 



6. Results in the model (solid lines, arbitrarily normalized) for topological cross sec- 
tions of charged particles with |?7| < 2.5 in pp collisions, compared with experimental 
data from ^0[. Upper curves and data correspond to ^/s = 0.9 TeV, those in the 
middle (multiplied by 0.1) to ^/s = 0.5 TeV, and lower curves and data (multiplied by 
0.01) to = 0.2 TeV. 

7. Yields per unity of rapidity at central rapidity, as a function of the number of 
wounded nucleons, for A, S~ and + Q'^ (left), and for p, A and S+ (right), for 
pPb collisions and four different centralities in PbPb collisions at SPS energies. Full 
lines represent our calculation with string fusion, and dashed lines with fusion and 
rescattering. Experimental data are from the WA97 Collaboration [Q. 

8. Results in the model (dotted line: without fusion, dashed line: with fusion, solid 
line: with fusion and rescattering) for strange baryon production in central PbPb 
collisions (5 % centrality) at SPS compared with experimental data from the WA97 
Collaboration |Q (triangles) and the NA49 Collaboration (squares). 

9. Results in the model (filled circles: with fusion, open circles: with fusion and 
rescattering) for the inverse exponential slope of the ttit distributions at midrapidity 
of different particles versus the mass of the particles in central (5 % centrality) PbPb 
collisions at SPS, compared with the experimental data of the WA97 Collaboration 



(3.5 % centrality, open squares). We also present our predictions for the same collisions 
at RHIC energy with fusion, filled triangles, and with fusion and rescattering, open 
triangles. 

10. Results in the model for the p — p rapidity distribution in central (5 % centrality) 
PbPb collisions at SPS (a), solid line), and RHIC and LHC (b)), dashed and dotted 
lines respectively), compared with experimental data at SPS |]72| . 

11. Results in the model (dotted line: without fusion, solid line: with fusion, dashed 
line: with fusion and rescattering) for the p rapidity distribution in central (5 % cen- 



trality) PbPb collisions at SPS, compared with experimental data [74 . 



12. Results in the model for the rapidity distribution of negative particles in central 
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(5 % centrality) PbPb collisions at = 17.3 GeV per nucleoli, without fusion (dotted 
line), with fusion (dashed line) and with fusion and rescattering (solid line), compared 
with data from NA49 |7|. 

13. Results of the code for the pseudorapidity distributions (upper plots), and the pr 
distributions for particles with \ri\ < 2.5 (lower plots), of charged particles in central 
(6 < 3.2 fm) AuAu collisions at y/s = 200 GeV per nucleon. In the plots on the left. 



solid lines are results with EKS98 ||33[ parametrization of parton densities inside nuclei. 



dashed lines with a parametrization as F2A |34], and dotted lines without modification 



of parton densities inside nuclei. In the plots on the right, solid lines are results without 
string fusion, dashed lines with string fusion, dotted lines are nucleon-nucleon results 
at the same energy, scaled by the number of wounded nucleons (344.6/2), and dashed- 
dotted lines are results with string fusion and rescattering. 

14. Results of the code for the pseudorapidity distribution (upper plot), and the pt 
distribution for particles with \ri\ < 2.5 (lower plot), of charged particles in minimum 
bias pp collisions at ^/s = 14000 GeV. Line convention is the same as in Fig. |^ (the 
dashed-dotted line is absent). 

15. Results of the code for the pseudorapidity distributions (upper plots), and the pt 
distributions for particles with \ri\ < 2.5 (lower plots), of charged particles in central 
{b < 3.2 fm) PbPb collisions at ^/s = 5500 GeV per nucleon. In the plots on the left. 



line convention is the same as in Fig. 13 left, but results have been obtained without 
rescattering and the dashed-dotted lines show results without semihard contribution. 
In the plots on the right, solid lines are results without string fusion, dashed lines with 
string fusion, dotted lines are nucleon-nucleon results at the same energy, scaled by the 
number of wounded nucleons (382.1/2), and dashed-dotted lines are results with string 
fusion and rescattering. 
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Tables: 





No fusion 


Fusion 


Ex])erinient 


charged 


7.89 


7.81 


7.69 ±0.06 


negatives 


2.95 


2.90 


2.85 ±0.03 


P 


1.18 


1.19 


1.34 ±0.15 


7r+ 


3.40 


3.33 


3.22 ±0.12 


7r~ 


2.69 


2.63 


2.62 ±0.06 


7r° 


3.73 


3.68 


3.34 ±0.24 


K+ 


0.31 


0.32 


0.28 ±0.06 


K- 


0.18 


0.18 


0.18 ±0.05 


A 


0.223 


0.231 


0.096 ±0.010 


A 


0.029 


0.033 


0.0136 ±0.0040 


P 


0.059 


0.070 


0.05 ±0.02 



Table 1: 
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No fusion Fusion Experiment 



p 


1.20 


1.21 


1.20 ±0.12 


7r+ 


4.04 


3.94 


4.10 ±0.26 


7r~ 


3.32 


3.23 


3.34 ±0.20 


7r° 


4.47 


4.38 


3.87 ±0.28 


X+ 


0.38 


0.38 


0.33 ±0.02 


K- 


0.25 


0.24 


0.22 ±0.01 


A 


0.245 


0.251 


0.13 ±0.01 


A 


0.045 


0.049 


0.020 ± 0.005 


P 


0.088 


0.100 


0.063 ± 0.002 



Table 2: 





No fusion 


Fusion 


Exi)erinicnt 


pS 


5.01 


4.86 


5.10 ±0.20 


pAr 


5.31 


5.12 


5.39 ±0.17 


pAg 


6.57 


6.28 


6.2 ±0.2 


pXe 


6.89 


6.56 


6.84 ±0.13 


pAu 


7.54 


7.16 


7.0 ±0.4 



Table 3: 
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NF F FR Experiment 



negatives 108.2 101.3 100.7 98 ± 3 



K+ 


9.7 


10.4 


10.8 


12.5 ±0.4 


K- 


7.1 


7.2 


7.4 


6.9 ±0.4 


A 


5.0 


5.9 


6.0 


9.4 ± 1.0 


A 


0.4 


1.1 


1.2 


2.2 ±0.4 


P 


0.82 


3.23 


2.80 






0.024 


0.186 


0.205 




2+ 


0.028 


0.097 


0.102 




Q- 


0.001 


0.007 


0.010 




n+ 


0.001 


0.005 


0.007 





Table 4: 





NF 


F 


FR 


Experiment 


A/A 


0.14 


0.41 


0.34 


0.128 ±0.012 




1.12 


0.52 


0.45 


0.266 ±0.028 


Q+/9- 


0.75 


0.88 


0.59 


0.16 ±0.15 


S-/A 


0.01 


0.06 


0.08 


0.093 ± 0.007 




0.08 


0.07 


0.10 


0.195 ±0.023 


Q/E 


0.05 


0.05 


0.11 


0.195 ±0.028 



Table 5: 
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RHIC(F) 


RHIC(FR) 


QCM 


Rafelski 


B-M 


LHC(F) 


LHC(FR) 


A/ A 


1.01 


0.90 


0.69 


0.49 ±0.15 


0.91 


1.00 


0.98 




0.96 


0.97 


0.83 


0.68 ±0.15 


1.0 


0.98 


0.95 


n+fn- 


1.00 


1.25 


1.0 


1.0 


1.0 


0.76 


1.03 


H-/A 


0.10 


0.15 




0.18 ±0.02 


0.13 


0.09 


0.25 




0.10 


0.16 




0.25 ±0.03 


0.14 


0.09 


0.24 


fl/E 


0.07 


0.26 




0.14 ±0.03 


0.20 


0.05 


0.40 


A/p 


0.40 


0.71 




2.4 ±0.3 


0.52 


0.35 


0.82 


P/p 


0.93 


0.90 


0.58 




0.90 


1.00 


1.04 



Table 6: 
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